Completely Quantized Collapse and Consequences 



3 



(N 



> 



Philip PearlCl 

Department of Physics, Hamilton College, Clinton, NY 13323 
(Dated: February 1, 2008) 

Promotion of quantum theory from a theory of measurement to a theory of reality requires an 
unambiguous specification of the ensemble of realizable states (and each state's probability of real- 
ization). Although not yet achieved within the framework of standard quantum theory, it has been 
achieved within the framework of the Continuous Spontaneous Localization (CSL) wave function 
collapse model. In CSL, a classical random field u>(x, t) interacts with quantum particles. The state 
vector corresponding to each tu(x, t) is a realizable state. In this paper, I consider a previously 
presented model, which is predictively equivalent to CSL. In this Completely Quantized Collapse 
(CQC) model, the classical random field is quantized. It is represented by the operator W(x, i) 
\f} , which satisfies [W(x, t), W(x', t')] — 0. The ensemble of realizable states is described by a single 

state vector, the "ensemble vector." Each superposed state which comprises the ensemble vector at 
time t is the direct product of an eigenstate of W^(x, t'), for all x and for < t' < t, and the CSL 
£N) i state corresponding to that eigenvalue. These states never interfere (they satisfy a superselection 

rule at any time), they only branch, so the ensemble vector may be considered to be, as Schrodinger 
put it, a "catalog" of the realizable states. In this context, many different interpretations (e.g., many 
worlds, environmental decoherence, consistent histories, modal interpretation) may be satisfactorily 
applied. Using this description, a long standing problem is resolved, where the energy comes from 
that particles gain due to the narrowing of their wave packets by the collapse mechanism. It is 
shown how to define energy of the random field and its energy of interaction with particles so that 
total energy is conserved for the ensemble of realizable states. As a byproduct, since the random 
field energy spectrum is unbounded, its canonical conjugate, a self-adjoint time operator, can be 
discussed. Finally, CSL is a phenomenological description, whose connection to, or derivation from, 
more conventional physics has not yet appeared. We suggest that, because CQC is fully quantized, 
it is a natural framework for replacement of the classical field ui(x, t) of CSL by a quantized physical 
\Q ' entity. Two illustrative examples are given. 
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I. INTRODUCTION 

Quantum theory as a theory of measurement was promulgated by Bohr[|, and this point of view has its adher- 
q-i ents todayH. But, the success of the theory suggests that it could be more, a theory of reality. Einstein wrote 
Schrodinger |3| that, if there were no more complete description of nature than that envisaged by Bohr, "...then 
. £h ! physics could only claim the interest of shopkeepers and engineers; the whole thing would be a wretched bungle." Bell 
made much the same point, writing^] " To restrict quantum mechanics to be exclusively about piddling laboratory 
Jh \ experiments is to betray the great enterprise." 

The lack in standard quantum theory has been called the "measurement problem," but I prefer to call it the 
"reality problem." I believe it is most precisely characterized as a failure to give a satisfactory response to the 
following operational requirement for a well-defined quantum theory of reality. You are given, at time t = 0, a state 
vector for an arbitrary physical system, together with its Hamiltonian, but no other information. For t > 0, specify 
(i.e., give a procedure for how to identify) the realizable states and their probabilities of realization. 

Not only Bohr's version of quantum theory, but more recent sophisticated attempts, such as many worlds 
environmental decoherence 0, consistent histories^], have not succeeded so far at providing such a specification^ 
They are successfully applied to many important and interesting examples, but these examples require additional 
information {this is the nature of the system, that is the apparatus, there is the environment...) which lie outside the 
theory. 

The reality problem has been solved by enlarging quantum theory so that it describes state vector collapse as a 
physical process. The Continuous Spontaneous Localization (CSL) model [ij, ITol Hlf provides a modified Schrodinger 
equation. It describes how, in addition to their usual interactions, the particles of the world interact with a posited 
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classical field w(k, t) resulting in a nonunitary state vector evolution (which changes the norm of the state vector). 
What are the realizable states? Each such state vector, corresponding to a different u>(x, t) (drawn from the class 
of white noise "functions), is realizable. What are the associated probabilities? CSL provides a second equation, 
the prescription that the probability of occurrence of a state vector (the probability of occurrence of its associated 
u>(x, t)) is proportional to its squared norm. In this straightforward manner, the operational requirement cited above 
is satisfied. Each of these CSL states accurately describes (except for occasional brief moments) the macroscopic world 
we see around us. The predictions agree with all presently performed experiments (but there are testable differences 
from standard quantum theory). The CSL model is reviewed in section II. 

I have previously presented a modeljl, 0, 0] called the "Index" model, but here called the "Completely Quan- 
tized Collapse " ( CQC) model. It is reviewed in section III. It utilizes a formalism invented to model continuous 
measurements 1 13j . CQC is completely equivalent to CSL in its specification of the realizable states and their proba- 
bilities, in the following manner. 

In CQC, the classical field is replaced by a quantized field W(x,t) which, like the classical field w(x,t) of CSL, 
commutes with itself, [W(x, i), W(x.' , t')] — 0. There is a conjugate field, II(x, <), which also commutes with itself, but 
which satisfies [W(x, t), II(x', t 1 )] ~ iS(t — t')S(x — x'). The basis of eigenstates of W(x,t) for every x, t, is denoted 
\w). It satisfies W(x.,t)\w) = w(x,t)\w), i.e., there is a different eigenstate for every possible u>(x, t) function. There 
is a single state vector, called here the "ensemble vector," which evolves unitarily (the Hamiltonian is given). When 
the ensemble vector is expanded as a Schmidt decomposition in the "preferred basis" \w), the particle state in the 
direct product with \w) is the CSL state associated with iu(x, t): schematically, 



The usual quantum rule for probabilities (the squared norm of a state in the superposition) gives the CSL probability 
of realization of that state. 

The measurement (reality) problem in standard quantum theory has often been phrased in terms of difficulties 
associated with the ill-defined collapse postulate, nowhere more succinctly than Bell's Boolean phrase "And is not 
Or," i.e., Schrodinger's equation gives a sum of terms but we observe one term or another. CSL provides a resolution 
of the problem by individually, dynamically, providing each term we can observe. However, CQC gives a sum of 
terms: how can that be considered as a resolution of the problem? 

A colleague once observed to me that construction of CQC goes against everything I have worked for, the resolution 
of the reality problem by dynamical collapse. However, what I have worked for is to make a well-defined quantum 
theory of reality i.e., one which satisfies the aforementioned operational requirement. CSL satisfies this require- 
ment, but so does CQC. The terms in the superposition which make up the ensemble vector are precisely defined (by 
the orthogonal preferred basis \w)), do not interfere (they obey a superselection rule), are macroscopically sensible 
(because the associated CSL states are) , and occur with the correct amplitudes (which give the correct probabilities). 
In his famous "cat paradox" paper |l5j. Schrodinger referred to the state vector as the "prediction catalog" or the 
"expectation catalog" because, in Bohr's quantum theory, it lists what can be expected to occur if a measurement is 
performed. The CQC state vector may, likewise, be called a "reality catalog." 

However, the justification here is greater than Schrodinger's. His usage requires judgment outside the theory as to 
when the state vector has evolved sufficiently to be considered a catalog, and as to which are the appropriate states 
in the superposition (i.e., which basis) to consider as making up the catalog. In CQC, these criteria lie inside the 
theory, since the ensemble- vector may always be considered a catalog, and the appropriate states are labeled by the 
preferred basis \w). Roughly speaking, CSL provides a vertical catalog (Or), and CQC provides a horizontal catalog 
(And) , but it is the same catalog. If one has a precise rule for identifying a superposition of states at time t which 
never interfere thereafter, it is a matter of indifference whether one removes them from the superposition and follows 
their separate evolutions that way, or leaves them in the superposition and follows their joint evolution that way. 

As previously pointed outy|, one may regard CQC as a model for a satisfactory resolution of other attempted 
interpretations of quantum theory^^. The superposed states making up the ensemble vector could be thought of as 
many worlds (although it would seem that there would then be no reason to insist upon such an interpretation). If the 
VF-field is regarded as the environment, then CQC could be regarded as describing environmental decoherence. The 
evolving sup erposed states making up the ensemble vector could be taken to represent a single family of consistent 
histories [lj]. The Schmidt decomposition of the ensemble vector is precisely what the modal interpretation requires. 
What is currently missing from these approaches, and what CQC provides, is a subsystem with a well-defined "uni- 
versal basis" (here the \w) basis vectors) whose Schmidt decomposition partners make macroscopic sense. One might 
hope that such a basis may be found arising in a natural way from a physical entity, whence perhaps CSL/CQC and 
all these interpretations may merge (see the last two paragraphs of this introduction for further remarks on this). 
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Now I turn to the topics which will occupy the rest of this paper. Because collapse narrows wave functions, the 
energy of particles in CSL increases with time. This growth of energy of particles provides for experimental tests of 
the model[ll| but, from a theoretical point of view, it is desirable to have a conserved energy. Recently, I showed how 
to define energy associated with the w-field in CSL so that the first moment of the total energy (particles plus w-field) 
is conserved [l9j. However, energy conservation requires all moments of the energy to be conserved. The formalism of 
CQC facilitates that demonstration (which can be converted to CSL language with less transparency). Because the 
ensemble vector evolution in CQC is a Hamiltonian evolution, the Hamiltonian is a conserved energy associated with 
the ensemble. 

The Hamiltonian does not commute with W(x, t) since it is composed of the W and n-field operators. Thus one 
cannot associate a conserved energy with an individual realizable state labeled by eigenvalues of W . Nonetheless, it 
is of interest to see how energy fares for the ensemble of possible states, for example, how the ensemble energy of the 
H^-field can go negative to compensate for the energy increase of the particles whose wavefunctions are narrowed by 
collapse. The discussion of energy conservation takes place in section IV. 

Another example of the utility of the CQC model is considered in section V. Because the VK-field energy has the 
whole real line for its spectrum, it is possible to construct its self-adjoint conjugate, a time-operator T, built solely 
out of the W and n-field operators. One may consider that time should be intimately related to the realization of 
physical states, since without such a realization there are no events and without events there is no time. This, and the 
fact that, like the Matterhorn, it is "there" (i.e., very few viable physical theories possess an unbounded Hamiltonian) 
motivate looking at T. 

X"s ensemble probability distribution is examined. T's variance diminishes with time. T's mean, roughly speaking, 
is a "center of time," i.e., the average time over all space, and over time from to t, weighted by the square of the 
particle number density. Thus, if the particle Hamiltonian vanishes so the only particle evolution is collapse, which 
does not change the ensemble spatial particle density, the mean of T is t/2. If, however, the particles have a high 
clumping rate (for example, due to gravitational interaction), the spatial average of the square of the particle number 
density is largest most recently, and the mean of T is closer to t. 

One of the motivations of this work is to provide a bridge from dynamical collapse models to other interpretations 
which seek a well-defined quantum theory based upon the standard quantum formalism alone (although they may 
be satisfied with just accounting for observation, not reality), rather than the altered formalism of collapse models. 
Now, what I consider to be a pressing problem associated with CSL is that it is phenomenological. (Some have said 
"ad hoc," but that is not appropriate terminology since ad hoc means "for this case only," whereas CSL applies to 
all cases, all physical situations: in fact, the label "ad hoc" more appropriately belongs to these other interpretations 
to the extent that they require additional information appropriate to each physical situation). It could be that CSL 
dynamics may arise from a larger structure which generates quantum theory itself, as has been suggested by Adler[20|. 
In this paper, we consider another possibility. 

CSL depends upon interaction of a posited field's interaction with particles, but no identification of that field with 
a physical entity has been made. A purpose in discussing CQC is that, since it is formulated completely quantum 
mechanically, that may make it easier to identify the universal W-field with a quantized physical entity which arises 
naturally in another context. In that case, CQC would appear as based upon the standard quantum formalism alone. 

However, it is not so easy to find a candidate for a universal fluctuating field with the desired properties. There 
is experimental evidencejlg that the coupling of W is to the mass-density of particles, which strongly suggests 
a gravitational connection, a proposal which has a long history [22, H3|. This sugggests looking for a mechanism 
whereby gravitational excitations could naturally act like the W^-field. 

By way of illustration, Section VI contains two models for W with gravitational overtones, neither to be taken too 
seriously. In one, W is constructed out of many scalar quantum fields (of Planck mass), each interacting only for a 
brief interval (the Planck time) with the particle's mass-density. The other utilizes aspects of a previous gravitational 
classical model 23] of W and crudely pays homage to spin- network models [2^. It assumes that space consists of 
Planck- volume cells, each containing a bound spin 1/2 entity with two possible energy states equal to ±Planck mass. 
A spin occasionally breaks free to briefly interact with a thermal bath, as well as gravitationally with the particle's 
mass-density, before it is bound once more. In this case, W(x, t) is identified with the free mean spin, in a small 
volume surrounding x, during the small time interval surrounding t,. These examples suggest that what is required 
for W is a quantity which has the freedom to vary independently over all space-time and for which the state vector 
describes its prior space-time values . This is unlike a conventional field, which can vary independently (in magnitude 
and time derivative) only on a single spacelike hypersurface, and whose present values are all that is described by the 
state vector. 
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II. CSL 

A. The Simplest Model 

The simplest CSL model, which will supply most of the illustrative calculations in this paper, describes collapse 
toward an eigenstate of the operator A. The state vector evolution, governed by a function w(t) of white noise class, 
is given by 

\ll>,t) S =T e -^'{^A + (4A)-V( t ')-2AA]2 }|( ^ (1) 

In Eq.(l), T is the time-ordering operator, Ha is the (usual) Hamiltonian of the system, \<p) = |^>,0) is the initial 
state vector, the subscript S denotes the Schrodinger-picture state vector, and the parameter A governs the collapse 
rate. Eq. (1) is the solution of the modified Schrodingcr equation: 

d\ip,t) s /dt = {-iH A - (4X)- 1 [w(t) - 2XA] 2 }\^,t) s . 

However, only the solution form (1) shall be utilized in this paper. 

CSL is completed by specifying the probability that a particular w, lying in the interval (w(t), w(t) + dw), is realized 
in nature: 

P t (w)Dw= s {ip,t\^,t) s Dw. (2) 

In Eq.(2), and for all integrations, one may discretize time, i.e., with tj = jdt, each w(tj) is considered to be an 
independent variable (— oo < w(tj) < oo), and Dw = n t j=l{dw(tj)/[dt/2ir\] 1 / 2 }, with the eventual limit dt — > 0. 
In the "collapse interaction picture," which is what is mostly employed in this paper, the state vector is 

\il>,t) = expiHAt\ip,t)s (3a) 

= Te- {iXrl Iodt'[w(t')-2\A(t')] 2 ^^ ( 3b ) 

where (3b) follows from Eq.(l), and A(t') = cxp(iHAt)Aex-p(— iH^t). From Eq.(3a), the probability rule, Eq.(2), is 
the same, except that the subscript S is dropped. 

The density matrix describing the ensemble in the interaction picture follows from Eqs.(2), (3b): 

Pit) = J PtM^ J^j^ 1 = J Dw\^t)(^t\ (4a) 
= Te-^/o dt'lAUt')-A R (t')f l(j>){(l)l (4b) 

where, in an expansion of the exponent in Eq.(4b), A^t') (An(t')) appears to the left (right) of \<fi)((f>\, and is 
time-ordered (reverse-time-ordered) by T. 

B. Proof of Collapse When H A = 

If Ha = (or if the collapse rate is so fast that, over the time interval (0,i), we can well approximate A(t) « A), 
then the state vector dynamically collapses, as t — > oo, to one of the cigcnstates of A. Let \(j>) = ^2 n a n \an}, where 
A\a n ) — a n \a n ). An intimation of the collapse behavior in this case may be seen from the density matrix (4b): 



p(t) =Y,ana* m e-W^- a ™y\a n ){a m \, 

since the off-diagonal elements vanish as t — > oo, and the diagonal element magnitudes are the ones obtained from 
\<j)) using the "collapse postulate" of standard quantum theory. However, this ensemble density matrix behavior does 
not prove collapse in the individual case, since many different ensembles can have the same density matrix. 
Here is a proof of collapse not previously given. Eqs.(3),(2) become respectively: 



\^t) = Y,*ne-^- 1 ti dt '^- 2Xa ^\a n ) (5a) 

(5b) 



\<P,t) = ^ a „ e -( 4A )" 1 /o*'[-(0~2Aa„]^ |an 
n 

P t { W ) = ^| an |2 e -(2A)- 1 / t <i*'[-(O-2Aa n ] 2 
n 
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It shall now be shown that the only w(t)'s of non-zero measure are those for which the asymptotic (t — > oo) 
time-average of w(t) is 2Xa n . To see this, consider the joint probability of w and a = (2Xt)^ 1 f Q dt'w(t'): 

P t {w,a)Dwda = 2Xtda5[ [ dt'w(t') - 2Xta]P t (w)Dw, (6) 

Jo 

where Pt(w)Dw is the probability of w and 2XtdaS is the conditional probability, given w, that u>'s time-average is 
2Xa. To find the probability distribution -Ft (a) of a, one must integrate (6) over Dw. To do so, make the replacement 
8{z) = (27r) -1 J dui expiwz, and then the resulting gaussian integrals over Dw may be performed. Next, the (gaussian) 
integral over uj can be performed, resulting in: 

e -2Xt(a-a„) 2 

P t (a)da = Kl 2 r~FTT7 da t=^2 l a «l^( a _ a n )da. (7) 

Thus, the asymptotic ensemble probability is completely accounted for by a taking on only the values a n , and the 
probability of an a n is \a n \ • 

It remains to show that (5a) asymptotically describes collapse. Let w(t) satisfy the constraint C a defined as 
limT^ 00 (2AT) _1 L dt'w(t') — a. Taking the asymptotic limit of (6) and then integrating over an infinitesimal range 
about a followed by integrating over w, it follows from (6), (7) that 

/pa+e /• 
Dw / daPoc(w, a) = / DwP oc (w). 
Ja-e Jc a 

Applying this result to (5b) yields 

C a 

Since each term in the integral here is non-negative, if a ^ b, for each w(t') obeying the constraint C a (i.e., for each w(t') 
which asymptotically equals 2Xa plus a fluctuating term which averages out to 0), then exp — (2A)" 1 J °° dt'[w(t r ) — 
2Xb} 2 = 0. 

This result can be applied to Eq.(5a), where the gaussians are of the same form (except 4A replaces 2A). Thus, for 
each such w(t) which obeys the constraint C a , with a ^ a„, the state vector asymptotically vanishes: that doesn't 
matter since, by Eq.(5b), the probability of its occurrence vanishes also. When w(t) obeys the constraint C Qm , all 
terms in (5a) but the mth asymptotically vanish, \ip, oo) — >^ \a m ): collapse occurs. 



C. Full-Blown CSL 



The CSL proposal for a serious collapse model, applicable to all nonrelativisitic systems, involves multiple A's, each 
with its own w. That is, collapse is toward the joint eigenstates of an operator A(x) for all x, and w(t) is replaced by 
the random field w(x, t). The choice of A is[9ll23|: 

>«»)-E^/*«-^ rW «. , (.)e.(.). <s) 

In Eq.(8), ^i(z) is the annihilation operator at location z for a particle of mass m, (mo is the proton mass), and a 
(and A) are parameters from the seminal GRW model[25j, with chosen values a w lCP 5 cm and A w 10 _16 scc _1 . Thus, 
since ^](z)^(z) is the particle number density of the ith particle type at z, ^4(x) is essentially proportional to the 
mass density inside a sphere of radius a centered at x. The CSL dynamical equation which replaces Eq.(3b) is: 

IV, t) = Te-^ 1 /o* rf*'dx^(x,t')-2AA(x, t ')] 2 |0) (9) 

(where ^4(x, t') = exp(iHAt)A(x) exp(— iiJ^t)), i.e., time is replaced by space-time. The probability rule is unchanged, 
Pt(w)Dw — (ip,t\ip,t)Dw, but space-time is discretized for integration, so Dw = Y[j k {dw(xk,tj)/[dtdx/2iTX} 1 / 2 }, 
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III. CQC 

A. Simplest Model 

Start by introducing creation and annihilation operators b'((j), b(cu) for — oo < uo < oo, which obey the usual 
commutation rules, 

[6(w), = 0, [fit (w), &*(«)'] =0, [&(«), &t(a/)] = J(w -£</)• 
Then, the operator W(f) (the quantized equivalent of the CSL w(t)) and its conjugate H(t) may be defined as 

/oo 
^[e-^fcH + e^VH], (10a) 
- OO 

n(t) = z/(8ttA) 1/2 / dw[-e-^*6(w) + e^ t 6 t (w)] (10b) 



which are readily shown to satisfy 

[w(t),w(t')] = o, [n(i),n(0] = o, [w(t),n(t')] = w(t-f) (11) 

The vanishing of the first two commutators in Eqs.(ll) result from cancellation of the positive frequency contribution 
in Eq.(10) by the negative frequency contribution. (Incidentally, a free tachyonic field also has equal positive and 
negative energy spectra, so it too can be quantized to have vanishing self-commutator |2^). 

The joint eigenstate of W(t) (H(t)) is denoted \w) (\tt)), where the label represents a particular function w(t) 
(7r(i)) of white noise class. It satisfies W(t)\w) = w(t)\w) (II(t)|7r) = 7r(t)|7r)). The range of t is (—00,00), and the 
states are normalized so that D'w\w)(w\ = 1 (f^° D'Tr\Tr)(n\ = 1). D' differs from the D employed for CSL 
in that it involves the product of dw(t)'s for the infinite range of t. To discretize, write \w) = YiJL-oo \ w {tj)i with 

(M*j)Kfe)) = S[w(tj) - w'(tj)W2ir\/dt. 
It follows from Eqs.(ll) that, for an arbitrary vector |$), 

{wmm = l6^) {wm - (12) 

Here, S/Sw(t) is the functional derivative. In time-discretized calculations, 5/5w(t) — > (dt)~ 1 d/dw(tj). 

The vacuum state |0) satisfies b(cu)\0) = as usual. Since, from Eqs.(10), W(t) + 2i\H(t) is an integral over b(u>), 
Eqs.(10), (12) imply 

(w\W(t) + 2iAn(i)|Q) = [w(t) + 2XS/8w(t)](w\0) = 
and a similar equation for (7r|0), whose solutions are 

HO) - e -(4A)- 1 .ffL*'- 2 (O j (?r | 0) = e-A/^dt'^W. (13) 

The normalization in Eqs.(13) ensure that / D'w\ (w\0) \ 2 = J D'n\ (tt|0) | 2 = (0|0) = 1. 

The W- field Hamiltonian, the energy operator which generates time-translations of W(t), H(t), is 

/•oo />oo 

H w = / duju:b\Lu)b{uj) = / dt'W(t')U(t'), (14) 



which follows from Eqs.(10). (By adding an infinitesimal imaginary part to u) in Eqs.(10), one can arrange that W(t), 

II(t) vanish at t — ► ±00, so one can freely integrate Eq.(14) by parts, to put the time derivative on II. One can also 
change the order of W(t') and Il(i') in Eq.(14), since / dt'[W(t'), U(t')} = 0.) 
The "ensemble vector" of CQC in the collapse interaction picture is defined as: 

|V>,t) = Te- 2iX ti-lt"Mt")m")\Q)\ ( f ) ) (l 5a ) 

= [ D'w\w)(w\Te' 2t ^o dt " A (nn(t")^^ (15b) 

' D'w\w)Te- 2X ti ^(0(J/h(0) e -W-' J". dt'^(t') )0) (15c) 

= /' D'w\w)Te- {ixrl /o dt'k(*')-2AA(i')] 2 | ( ^) e -(4A)- 1 {r oo if i» 2 (t')+/ t ~ dt'» 2 (t')} ^ 15d ^ 
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A way to look at this is as a continuous measurement sequence. The change of the wave function at time t over the 
interval dt, according to (15a), resides solely in 

e -2i\A(t)i- 1 d/dw(t) e -(4\)- 1 dtw 2 (t) _ g-(4A) _1 dt[w(t)-2\A(t)] 2 

which has the form of the von Neumann measurement of A(t) by the pointer w(t). However, here the pointer 
uncertainty ~ (dt)" 1 is quite large. It is only the accumulation of a large number of measurements which can 
distinguish the eigenvalues of A(t) . One can visualize the collection of pointers for all t (and for all (x, t) in the case 
of full-blown CSL), each waiting, patiently, for its brief opportunity to perform its measurement. 

B. Equivalence of CSL and CQC 

The time-ordered term in Eq.(15d) is precisely the CSL expression (3b) for the state vector which evolves under 
the classical noise w(t') for < t' < t. The remaining exponentials are unaltered parts of the vacuum state, one for 
— co < t' < which never changes, and one for t < t' < oo which evolves. 

In standard quantum theory, a superselection rule [23] obtains at time t if thereafter the state vector may be written 
as a superposition of orthogonal states, each evolving independently of all the others (i.e., the matrix elements of 
the Hamiltonian between these states vanishes). In such a case, the state vector may be regarded as a "catalog," as 
discussed in the Introduction, merely a convenient device for listing the different states and their probabilities, since 
any state's evolution after time t can be separately described. 

In CQC, the ensemble vector admits a superselection rule at any time. Consider the superselected states at time 
t. The labels which distinguish them are the values of w(t') for < t' < t. The differently labeled states evolve 
independently thereafter. Each state corresponds to a different CSL state. The values of w(t') for -co < t' < 
and t < t' < oo are irrelevant as far as the particle physics over the interval (0, t) is concerned. The former never 
becomes relevant. For the latter, as time increases, say from t to t\, the values of w(t') for t < t' < t\ become 
relevant: a superselected state at time t splits into many superselected states over the interval (t,t±) (in a way that 
a many- worlds advocate might admire). Thus, as far as the particle physics is concerned, the ensemble vector may 
be regarded as merely a convenient device for listing the different CSL states, and this is the sense in which CSL and 
CQC are equivalent. 

The evolution (15) is unitary and guarantees conservation of probability if the usual quantum rule is employed, that 
the probability of realization of the state ~ \w)(w\tp, t) is D'w\(w\tfj,t)\ 2 . As far as the particle physics is concerned 
at time t, one may integrate the probability of a superselected set over w(t') for t' outside the range < t' < t: the 
result is the CSL probability P t (w)Dw. 

C. Full-Blown CQC 

Full-blown CQC corresponding to full-blown CSL described in Eqs.(8),(9), simply replaces time integrals by space- 
time integrals, e.g., 

W(x) = (A) 1/2 /(2tt) 2 f d 4 k[e- ihx b(k) + e ik - x tf(k)] 



|V, t) = Te- 2lX /o d 4 x'A(x')n(x') | j |^ 

replace Eqs.(lOa), (15a), etc. 

D. Nonwhite Noise CQC 

A CSL collapse evolution which generalizes Eq.(3b)|28| is 

|V, t) = Te" (4Arl dt'dt"[w(t')-2\A(t')]g(t'-t")[w{t")-2\A(t")] 

where g(t) is real, symmetric in t and has a positive spectrum, i.e., g(t) = J duj exp(iwt)/i 2 (w) with h{w) — h(—u>) = 
h*(ui). 
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A CQC counterpart requires denning W(t) and H(t) as in Eqs.(lO), except that the square bracket in the integral of 
Eq.(lOa) (Eq.(lOb)) is multiplied by h^ 1 ^) (h(u>)). This replaces the equation before (13) and Eq.(13a) respectively 
by 

[ / dt'g(t - t')w(t') + 2XS/5w(t)] (w\0) = 0, (w\0) = e - (4A)_1 dt'dt"w(t') g (t'-t")w(t") 



With the evolution equation the same as Eq.(15a), Eq.(15d) is replaced by 

lip, t) = Te" (4A)_1 dt ' dt "l w ( t ')^ A (t')ea, t (t')]g(t'-t")[ w (t")^2\A(t")e„, t (t")] ^ 

where @o,t(t') = 1 for < t' < t and = otherwise. For general H(lu), this equation is not the same as its CSL 
counterpart written above, except in the S-matrix limit where the interval (0, t) is replaced by (—00,00). This 
generalization shall not be considered any further in this paper. 



IV. ENERGY 
A. Pictures 

The CSL Schrodinger picture state vector (1) has a corresponding CQC ensemble vector: 

\i>,t) s = e- iHAt \i>,t) (16a) 
= Te- i i'o*'[^+2AAn(t')]| )|^ ) (16b) 

where (16b) follows from putting (15a) into (16a). Although A's time dependence in (15a) is removed in (16b), this 
is still not the Schrodinger picture for CQC's ensemble vector (nonetheless, we'll keep the S subscript on it), since 
the IT operator still has time dependence. It can be removed in what we'll call the sub-Schr6dingcr picture: 

\1>,t) SS = e- iH ^,t) s (17a) 

= e -«[H„+ir A +2A>UI(0)]| 0> |^ } (17b) 

Thus, since this is a unitary evolution, there is a conserved energy, the Hamiltonian in Eq.(17b). In the collapse- 
interaction picture it is 

H(t) = S Ha+Hw ^[H a + H w + 2XAU(0)}e-^ HA+Hw '> t =H A + H W + 2XA(t)U(t). 

It should be emphasized that, since H w and the interaction term depend upon II which has nonvanishing matrix 
elements between different cigenstates of W, the conserved energy belongs to the ensemble of realizable states, not 
to an individual state. A single collapse to the realized state generally will not preserve the energy distribution. (Of 
course, this is also the case in standard quantum theory with the collapse postulate.) Nonetheless, it is interesting to 
consider the probability distributions of the ensemble total energy and its parts. A conserved quantity gives insight 
to the behavior of the ensemble of possibilities. 

For example, one may see how the the ensemble average increase of particle energy, due to the localization mechanism 
of CSL, is compensated by the ensemble average decrease of the random field energy. And, laboratory experiments 
often do involve an ensemble of (approximately) identical initial states, for which therefore energy is (approximately) 
conserved for the ensemble of outcome states. Also interesting are circumstances where one can speak of conservation 
of energy associated with the ensemble of states that collapses to a distinctly identifiable outcome. For example, when 
considering a situation when the initial state vector is a superposition of particles in widely separated locations, the 
ongoing orthogonality of the differently located particle states ensures energy conservation for the ensemble of states 
which describe collapse to particles in one of the locations. 



B. Simplest Model, H A = 



We display here generating functions and probability distributions for the total energy and its components, for the 
model of section Ha, in the simplest case, H A = 0. 
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Using the expression for \ip,t) given in Eq.(15a), the generating function for the moments of the total energy is 
(?M|e-^ ( *^|V,i) 



E 

n 


a„ 


2^Q| e 2iAa„ J* dt'Tl(t')] e -i[H w +2\a n 


n(t)]^ e -2 J Aa„/ o t dt'n(t')]| ^ 


(18a) 


E 


a n 


2^Q|g — i[ff t „+2iAa„ J * dt'iII(i')+2Aa T 


>n(t)]/3 


(18b) 


n 

E 

n 


a n 






(18c) 


E 

n 


«n 


2(o| e -*^ e - 2jAa ™ /tf dt ' n (*') |0) 




(18d) 


E 

n 






9 dtV(t') 


(18e) 


E 


a™ 


2 e -(A/2)a=|/3|_ 




(18f) 



(18b) is the result of performing the unitary transformation on e~ iH W in (18a), the integration in the exponent of 
(18b) results in (18c) which confirms that the generating function is time- independent, (18d) expresses the unitary 
operator in (18c) in a different form, (18e) follows from H W \Q) = and insertion of the identity 1 = J D'7r|7r)(7r|, (18f) 
follows by discretization and performance of the gaussian integrals. 

The moments of the energy follow from Eq.(18f) using the formal Taylor series expansion, |/3| = + /3e(0) + 
(/3 2 /2!)2<5(0) + (/? 3 /3!)25'(0) + ... (where e(x) = x/\x\ for x ^ and e(0) = 0). Thus the mean energy is and the 
higher moments of the energy are infinite. The associated probability distribution more transparently entails this 
result: 



V(E) 



1 

2^ 



iE/3 



{i/>,t\e 



liM> 



= ^Ek 



Aa 2 /2 



E 2 + {\al/2f 



(19a) 
(19b) 



Although the total energy probability distribution is time-independent, the energy H w is not. The generating 
function here is 



(TMIe-^l^t) = ^| a „|2(0| e 2^a„/ o t d t 'n(t') e -2 l Aa„/ o t d t 'n (t '-/3)| 0) 



J2\ a n\ 2 J D'ne 



-2X f^ 7r 2 (t') 2iXa n [J* dt-r_-/]dt'K(t>) 



= K| 2 {[©(/? ~t) + 0(-t - /3)]e- Aa «' 

n 

+e(t-/3)e(t + / 3) e - Aa "i /3 i}. 



(20a) 
(20b) 

(20c) 



where (20a) utilizes (15a) and the time-translation capability of H w , (20b) involves going to the \ir) representation, 
and (20c) follows from performing the gaussian integrals. From the Taylor series expansion of (20c), one sees that the 
mean value of H w is zero and the higher moments are infinite. 

The probability distribution following from Eq.(20c) confirms this result: 



V W {E) 



(il>,t\e 



-Xatt 



-iH w 



|lM> 



ttS(E) 



sin Et E sin Et — Aa 2 cos Et 
+ 



E 



+ 



E 2 + (Aa 2 ) 2 
Aa 2 



E 2 + (A 



<) 2 \ 



(21a) 



(21b) 



Eq. (21) shows how the VF-field energy starting at t = with the vacuum energy E = 0, evolves toward a stationary 
energy distribution associated with the completely collapsed state. (Note that the asymptotic value of (21b) is not 
the same as (19b), although it is similar in form.) 
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The interaction energy's probability distribution is 



1 r°° 

Vj(E) = — / e iE ^{^t\e- iB ^,t) (22a) 
= JL r e lE ?Y\a n \ 2 [ DVe-^oo^^e- 2 ^"-^ (22b) 

27r J-oo „ 7 



-OO 

2 



= \- Fn| p -gV(2Aa'„/dt) (22c) 

V y/2n\al/dt 

Since the interaction energy is proportional to the white noise operator II(i), it has a white noise energy spectrum, 
with all energies equally likely. 

C. Simplest Model, H A 

1. Moment Generating Functions 

As in the previous section, also in this more general case, the random held variables can be integrated out of the 
expressions for the moment generating functions. For the total energy, corresponding to Eq.(18), one obtains: 

(i>,t\e- iH W\i),t) = M(Q| e -*[ H A+ff.+2An(o)A]/J| O )|0 (23a) 

\ e -i[H A +H w ]0 Te -2iXJidt'n(t')A(t')^^ (23b) 

|e -iH A /3 /' £) ' 7re -2A/- oo dt' 7 r 2 ( t ') Te -2iA/ ( fd i V(i')A(t')|0) ( 23 c) 

= (0| e - i ^/» Te -(A/a)/ o "«lfA»( t '}|0) ( 23d ) 

= (0| e -^+( A / 2 ) A2 ( )]^). (23e) 

((23a) is expressed in the manifestly time-independent sub-Schrodinger picture, and the remaining steps are ones 
taken before.) For Ha, the moment generating function is 

(iP,t\e- iHA(3 \ij,t) = (^|( |Te 2iA i'o*' n (*')^i.(t') e -w A ^ e -2iA/ t dt'n( t ')^(t')|o)| ( / ) ) (24a) 

= {^T{e-^ 2 y^ dt ^ A ^- A ^\- lH ^}\<t)) 1 (24b) 

where Ar (Al) acts to the right (left) of exp— iHaP, and is time-ordered (time-reverse-ordered) by T. 
For H w , the calculation is 

(V>,t|e- iff «'V,t) = (^|(0|Te wx ^««' n (OAi i (*') e -w,/J e -2iA/*«a'n(f)>i Jl (f)| ^^ (25a) 

= (^KOl^e 2 ^^" d *'n(*'+/3/2) J 4 i (i') e -2TA/ t dt'n(t'-/3/2)^ J7 (t') | ) |^,) (25b) 



Te 2iA j;^ 72 dt'7r(t')^i(t'-/3/2) e -2iA f*~% 2 dt'7r(t')A R (t'+/3/2) |^ 

= (^ire-^/^/o'dt'^icto+^ct')] [ (/3 _ *) + e(-/3 - 1) 



2. First Moments 



Here is shown conservation of energy in the mean, reproducing more simply a previously obtained result (all that 
that had been proved about conservation of energy for CSL19]). Equating the terms proportional to the first power 
of j3 in Eq,(24b) produces 
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To find the first moment of H w , the exponent in (25d) must be written for small (3: 

rt-\P\/2 pt 

/ dt'A L (t'-(3/2)A R (t' + (3/2) = / dt'A L (t')A R (t') 

J\P\/2 JO 

-(\/3\/2)[A L (t)A R (t) + A L (0)A R (0)] + (/3/2) f dt'[A L (t')A R (t') - A L (t')A R (t')} + ... (27) 

Jo 

The zeroth and first order terms in the Taylor series expansion of (27) are the same as (27), except that the second 
term vanishes, since d\(3\/dt\p=o = e(0) = 0. The step-functions in (25d) likewise do not contribute to the Taylor 
series expansion's first two terms, so we get from (25d): 

(V,WV,*> = i (A/2)(0|r e -( A / 2 )/o t ^'[^(*')-AH(t')] 2 

■/ dt'[A L (t')A R (t') - A L (t')A R (t')]\<i>) (28a) 
Jo 

= (\/2)(cf>\T [ rft'e-CVa)^' d*"[^(*")-A«(*")] a [A(t'), i^4(t')]|^) (28b) 
Jo 

= {\/2)(<t>\T I dt'e-^ fo dt "^ t "^- AR{t "^[A(t') 1 [A(t') 1 H A ]]\<j } ) (28c) 
Jo 

In going from (28a) to (28b) we have utilized 

Te -(V2)/ t ^[A i (t')-A«(t')] 2 SL( i l)Cij( i l)=Te -(A/2)/ ^ dt'[A L (t')-A R (t')f B{fi)C{jti) 

for t\ <t and arbitrary operators -Bl(£i), C R (t\) (the subscripts R, L may be dispensed with when operators having 
the same time argument are adjacent and in the correct order). 

Since the first moment of the interaction energy vanishes, from Eqs.(26), (28c) follows conservation of energy in the 
mean: 

(d/dt)(i>,t\H w \^,t) = -(d/dt)(i),t\H A \i,,t) 

= (A/2)(0|Te-( A / 2 ^o dt'[A L (t')-A R (t')f [ A {t)^ HA ]]^) ( 29 ) 

3. Large t 

The large t-behavior of the probability density distribution of H w is due to the third term in the bracket of Eq. (25d) . 
Because this term's exponent is largest for small (3, and also because expiE(3 oscillates rapidly for large /3, the small 
(3 approximation of (25d) should be a good approximation to the large t limit. Putting (27) into (25d) results in 



1 f°° 

V W {E) = — e iE P(ij,t\e- iH ^\i>,t) (30a) 

2n J-oo 

— / d(3e iEfi '((j)\Te- {X/2) /o d *'^^(*')->lK(t')] 2 e -(A/2)|/3|[A 2 (t)+A 2 (0)] e -i(A/2)/3/ 'dt'C I , K (t')|^ ( 30b ) 
27T J_ t 



= U\Te-^Iodt'[A L (t')-A R (t')fi \[Ai(t)+A\Q)]/2 

*[E- (A/2) /„' dt'C LR (t')f + [\(A*(t) + A^))/2Y 



(30c) 



where C LR (t') = A L (t')[A R {t'),H A ] - [A L (t'),H A ]A R (t'). Note that, when H A = 0, the large t-limit of (30c) is the 
same as (21b). Also, note that, when H A ^ 0, the mean value of E given by (30c) is the same as (28c). 

For example, in the case of a free particle of mass m which collapses toward a position cigenstate, where H — p 2 /2m 
and A = x, the term in (30c) [E - (A/2) / Q dt'C LR {t')} 2 = [E + Xht/2m] 2 . Because the ensemble energy is conserved, 
the mean VF-field energy decreases with time to compensate for the increase of the particle's mean energy = \ht/2m 
due to the narrowing of wavepackets by the collapse mechanism. 
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V. TIME OPERATOR 

Unlike the usual Hamiltonian which is semi-bounded[29(, H w has an unbounded spectrum, providing an opportunity 
to define a self-adjoint time operator: 



T = 



/£*'[(4A)-^ 2 (t') + An2(i')] ■ 



(31a) 
(31b) 



where [B,N] = 0. It is simple to show that [H W ,T] = i, where H w is given by Eq.(14). The action of T on |0), is 
undefined, but one can take T = B/(N + e) so that T|0) = 0/e = 0. Of course, T plus any function of H w and N is 
also conjugate to H w , but (31) is the simplest choice. 

A. Eigenfunctions 

T"s eigenvalues are highly degenerate. A complete set of (unnormalized) eigenfunctions in the \w) basis is: 

/oo 
dwe^(u;)] n ^\0) ~ JJ[ic(i)-2AVH n(t) e" (4Ar ( ° 

t 

where n(t) = 0,1,2... and the H n are Hermite polynomials. The associated eigenvalues Y^t^-oc t n (t)/ Ylt=—oo n W 
have a "center of time" form. 

B. Moment Generating Function 

To find the moment generating function for T, first observe, using (15a) and (10b), that 



-i/3T 



T 



1 + f e - Wn (V2^) n/2 r f dt'A R {t') r dwe^'tfiu)] 71 



|0) e -(V2)/o df,A H(*')| 



Since exp—(if3B/n)b^(uj)exp(i^B/n) = tf(u + 0/n), one gets: 



(A/27r)"/ 2 r [ l dt , AL{tl) r due-^'biio)] " 



i + y (A/27 ; )K/2 r /' dt'^(i') / ^(—""^^fttfw)] n 



|O>|0> 



ij-g-CA^/o'di'IA^CtO+A^t')] 



(A)' 



(32a) 
(32b) 



As a check, when /3 = 0, (32b) becomes: 

(0|re~ (A/2) ^ d *' [A « (t ' )+A ' (t ' )1 e A ^ dt '[^( t ')^i.(*')]|^ = ^|7- e -(A/2)/ f dt'[A fl (t')-A L (t')] 2 |^ = 1< 
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C. First Moments 

Now, extract the first two moments of T from (32b). From the coefficient of /3 in (32b) we obtain: 

f*dt't'A L (t')A R (t') 
f*dt'A L (t>)A R (t') 

^ e -(\/2)p dt'[A L (t')-A R (t')] 2 _ e -(X/2) f t Q dt'[Al(t')+Al(t')]^ |^ 

To get a feel for the meaning of (33), consider the case Ha = 0: 



l- fe | 2 4^[l-e- A ^] = |[l-E 



a k \ 2 e~ xt < 



(33) 



(34) 



Thus, asymptotically, (ip, t\T\ip, t) approaches the "center of time," a weighted average (here uniformly weighted) of 
the time over the interval (0,t). For Ha ^ much the same behavior prevails, according to Eq.(33). For example, 
set A R (t') = A L {t') = At' s . Then, from (33), asymptotically, (tp,t\T\tp,t) -> t(2s + l)/(2s + 2): for s large, the mean 
value of T is almost t. 

For full-blown CQC, the expression that replaces (33) is the same expression with A(t') — ► A(x, i') and dt' — ► dxdt'. 
A (given by Eq. (8)) is w the mass density in a sphere of radius a. Then the appropriate "center of time " expression 
appearing in the equivalent of (33) is 



f dt'dxt'A L (x,t')A R (x,t')/ [ dt'dxA L (x,t')A R (x,t'), 
Jo Jo 



i.e., t', roughly speaking, is weighted by the squared mass density integrated over all space and over time from to 
t. Should there be a high rate of agglomeration of matter, e.g., due to gravitational clumping, then (ip,t\T\ip,t) will 
approach t. It should be emphasized that gravitational collapse has this effect, not wave function collapse, as Eq. 
(34) attests. 

The second moment of T is likewise found from (32b) : 



<<MT 2 |V,t) = 



IT 



J*dt't'A L (t')A R (t') 



£dt'A L (t')A R (t') 

.[ e -(\/2)ftdt'[A L (t')-A R (t')] 2 _ e -(A/2) /„' dt'[Al(t')+A 2 R (t')}^ + f(Z LR )]] 



+ 



f*dt't' 2 A L (t')A R (t') (A/2) ( t dt > lA l(t')+Al(t')] f(z )) 
f*dt>A L (t')A R (t>) JK 7 



(35) 



where Z LR = A f* dt'A L (t')A R (t') and f(z) = f* dz'[cxpz' - l]/z'. For the example A R (t') = A L (t') = At' s (H A = 
corresponds to s = 0), Eq.(35) asymptotically becomes (using f(z) — ► z~ 1 expz): 



(^t\T 2 \^,t) 



t(2s + 1) 
2s + 2 



1 + 



1 (2s + 2) 2 ^ K' 2 



Ai 2s+1 2s + 3 — a 

k 



(36) 



so the fractional deviation decreases over time, AT 2 /T 2 <~ i (2-S+ 1 ), 



D. Large t 



More generally, an approximation to the large t probability distribution of T may be had by employing the small 
(3 approximation in (32b), 

[( dt'A L {t')A R {t')e- m '' n ] n « [[ dt'A L (t')A R (t')(l-0t'/n)] n 
Jo Jo 

r t 

w [ / rft'A I ,(t')A iJ (i / )]"e" i/3 -''o d *'*' A i(*')^K(*')//o ^'^(tO^fiC*')^ 
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with the result 



V(t) = (2tt)- 



(<j>\Te- (x/2)I ° dt'[A L (t')-A R (t')] 2 s 



f*dt't'A L (t')A R (ti) 
J*dt'A L (t')A R (t') 



(37) 



To summarize, there is a natural clock associated with the random field, whose time reading depends upon global 
dynamics and whose deviance decreases with time. In future work, rather than collapse depending upon time, one 
may look to express the attractive idea that the notion of time may depend upon collapse, since no collapse implies 
no events implies no time. 



VI. TWO EXAMPLES 



Collapse models are phenomenological, a guess (based upon a perceived flaw in standard quantum theory) that 
nature utilizes this description of dynamical behavior. To be generally accepted, such models need experimental 
verification j^Ol and grounding within a larger physical framework. One of the main aims of this paper has been to 
emphasize that CSL can be formulated as CQC, within an essentially quantum framework, in order to encourage 
thoughts about how collapse might arise from a physical structure which has a natural quantum theory expression. In 
that spirit, this section contains two examples of CQC dynamics arising from other quantum conceptual structures. 



A. Quantum Fields 



Consider a collection of scalar quantum fields, associated with particles of large mass M which, for definiteness, 
take to be the Planck mass, so large that Vk 2 + M 2 « M for any so-far experimentally achievable k. Each field (more 
precisely, its time derivative, which is its conjugate field) interacts with the mass density A(x, t) (given in Eq.(8) et. 
seq.) over a brief interval which, for definiteness, we shall imagine is the Planck time r = M _1 . Call the field ^((x) 
(and call its conjugate II t (x)) which is responsible for the interaction over the interval centered at time t\ 

d> t (x) ee C — — \-j= J dk{e^~ M % t (k) + e ^- M %l(k)} (38a) 

n t (x) ee C-'j^j-^ J dk[-e^- M % t (k) + e-^- M %\(k)} (38b) 

where the constant C is yet to be chosen. 

Define the eigenstates by &t(x-)\<ft) = <Pt(x-)\<ft) and the associated vacuum state by 6 t (k)|0) t = 0. Then, Eqs.(38) 
imply 

[^(x) + ^)K^( X )I°>* = °- M*)|0>* = e-^ dx ^ x >\ (39) 

The joint vacuum state for all the fields is |0) ee Y[ t |0)t and a joint eigenstate of all the fields is \(p) = Y[t I'/'tC*))*- 
To make their scalar product close to the scalar product in Eq.(13) (suitably generalized to space-time), choose 
C = MV2X so, from (39), 

The evolution equation is chosen as: 

{<p\i/>,t) = Te- 24A n< = o/^'nWx)A(x', t ')| O )|0) we -2iA/ t t , =o dt'/dx'A/n t ,(x')A(x',t')| 0) |^ (41a) 

= Te" (4Arl ^'=o T / dx'[ w (x')-2AA(x' \t')] 2 |0) e -(4A)- 1 E° = _ =o +i:^ =t ]^/d(x')^(x') ( 41b ) 

To complete the demonstration that this model is equivalent to CQC, identify the fields 

W(x, t) = $ 4 (x), ft(x, t) ee Mn ( (x) (42) 
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from which follows 

[W(x, t), n(x', t')] = 5(x - x.')5 t ,t>/T w (5(x - x ')<y(t - *')• 

It is immediately seen, using Eq.(42) for n(x, f), that the statevector evolution equation (41a) is identical to Eq.(15a) 
(generalized to space-time). 

To see that the expressions for W and II in terms of particle creation and annihilation operators closely approximate 
the similar expressions for W and II, note that 

[e- lMt 6 t (k)/V7, e lMt b\ (k)A/r| - 5(k - k')5 t<t ,/T « S(k - k')S{t - t'). (43) 

Therefore, define 

b{k,uj) = (2Tr)- 1 ' 2 J dt'e^'e-^'bvik!)/^?, (44) 

so [6(k, ui), b'(k', u)'] — S(k—k')S(u> — u)') follows from Eqs. (43), (44). Inserting the inverse of Eq. (44), for b t (k) in terms 
of b(k 7 lu), into Eqs. (38), (42) gives for W and n expressions whose form is identical to (10)'s expressions (generalized 
to space-time) for W and II. 



B. Many Spins 

The model described here borrows from a gravitationally based semi-classical collapse mo del [23^ and spin-network 
models of space-time [24j , without being faithful to either. Suppose that space consists of "Planck cells" (each of 
volume £ 3 , where £ is the Planck length), each containing a spin which is normally rigidly enmeshed with the other 
spins, but which occasionally breaks loose. Then, for a duration equal to the Planck time r, it interacts both with 
a thermal bath, (characterized by the energy /3 _1 ) and gravitationally with nearby particles (more on this shortly). 
Following this brief dynamics it is frozen, no longer interacting with either bath or particles, once more enmeshed 
with other spins. 

Characterize the spin in cell i located at x^, which interacts over the interval (t — r,t) and then freezes at time 
t, by <Tj ;t (a Pauli z-spin matrix) and mass /icr i t , where p is the Planck mass. Suppose that a spin- mass at x^ 
interacts gravitationally only with the particles in a sphere of radius a around it, and that the interaction is as if the 
particles are uniformly smeared over that sphere, with resulting mass density represented by the operator /o(xi). For 
simplicity, take the particle Hamiltonian to vanish, so that p is independent of t. Then, since any initial statevector 
for the particles can be written as the sum of joint eigenstates of p{x.i) for all x^, calculations for the statevector 
need consider only the special case of one such term, for which p(~x.i) is a c-number, with the more general case being 
reconstructed afterward as the sum. Then, for one such term, the Hamiltonian for the gravitational interaction of 
all the spins which have interacted with the particles over the interval (0, t) (up to a numerical factor which, for 
simplicity, is taken equal to 1) is 

t 

H = -Gfia 2 o-i,t'P(xi)- (45) 

i,t'—r 

Because all the spins act independently, focus upon a group of spins excited over the time interval (t',t' + At) 
which lie within a volume AV. Take At and AV to be large {At/r >> 1, AV/£ 3 >> 1), large enough to contain TV 
interacting spins where N >> 1, but still small on the scale of particle physics (e.g., p(xi) is essentially constant over 
AV^). Denote by p the probability that a spin in a cell is freed up to interact over a time interval r so, on average, 

(46, 

For simplicity, we shall take (46) to be the exact expression for N, not just for its average. 

Now consider the (interaction picture) statevector \ip, t). It may be written as the direct product of states associated 
with all volumes AV and all time intervals At within (0,i). The state which is the contribution of one such volume 
and interval has the form 

l l 

M = ■- Yj \si-SN)\hn,th, Sl ...s N )C Sl ... SN (p). (47) 
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The sum in Eq.(47) contains all the spin states entangled with the associated bath states which interacted with the 
spins and brought them to thermal equilibrium: assume for simplicity that the bath states are mutually orthogonal. 
Then, |x)'s contribution to the density matrix describing the spins alone is 



Tr bathlxKxl = H ■- H \ s i-- s n)(si--s n \\C Si ... Sn {p)\' 



Si=— 1 s. 



On the other hand, because these states are all in thermal equilibrium, using Eq.(45)'s Hamiltonian for these spins, 
the thermal density matrix is: 

■•• \si-sn){s 1 ...s n \— 1 -j ^— . 

si — — 1 s^ — — l Z-^/si — — 1 " ' Z-^sjy — — 1 

Thus, by comparison of these two equations, 

\ C ^~°M 2 = w w jG^JTZ^ ' (48) 

Eventually, W associated with this space-time region will be taken ~ S = X)i=i a i- The statevector which is the 
sum of all states in (47) which have the same eigenvalue of S and which is normalized to 1 is: 

C Sl ... SN (p) 



\s) = ^ \si-..s N ) |bath,si. ..s N ) 



{ Si },Y: si= s \zJ2{ Sz },y: s i=s \c Sl ... SN 12 

so, using (48), 



(s\x) 



-|l/2 P f3G Pt ia 2 Y;? =1 s,-\ 1/2 



(49) 



Evaluation of (49) is well known but, for completeness, it is sketched here. The partition function in the denominator 
of (49) is (writing C — Gppa 2 and taking N even) 

Tr^ cs =y e ^( N - 2k ^UN-k)l = (2 cosh (3C) N , (50) 
^— ' k 

k=0 

so the expression in (49) under the square root is 

t x? _ |Csi - Sn(p)|2 = (2 cosh /3C) W p - s)/2] ! p + a) /2] ! (51&) 

1 g-Is-ATtanh^ClVaWcosh 2 /3C (51b) 

V27r7Vcosh 2 /JC 

' , e -[s-N0Cf/2N^ (51c) 



V2^V 

In (51b) the well-known gaussian approximation to the binomial distribution has been employed[3l] which is valid for 
large N and for s within a few standard deviations of its mean. The approximation in (51c) depends upon validity of 
the inequality 

f3G Pf ia 2 « 1 (52) 

which shall need to be checked. 

Now it is possible to put together the contributions of all volumes AV and all time intervals At within (0, i) to 
construct the statevector \ip,t). Define |s) = JIava* \ s )i,ji the normalized basis vector which is the joint eigenvector 
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of S(x.i,tj) for all Xj and for < t j < t. The scalar product (s\ip,t) is the direct product of the expressions (49) so, 
from (51c) (without the normalization factors, which are tucked into the element of integration): 

<J|*,t) = e^O-^k^-J**-""**? (53a) 

e -4~ 1 N(0G t J.a 2 ) 2 -Ei i , tj = [si,j/N(3G t j.a 2 -p i ] 2 
~ g-r'W^TX/JG^) 2 / dx/ ' «'[ S '(x,t')-p(x)] S ! ( 53c ) 

where (53c) follows from (46), and we have defined: 

s'(x,i) = Si.j/NpGfia 2 . (54) 

Now, compare (53c) with CQC's comparable expression. Note from (8) that A(x) — m ~ 1 a 3 / 2 /9(x) (up to a numerical 
factor which shall be ignored) so the CQC expression comparable to Eq.(9) is 

(w\lp,t) =e -(^r 1 I^J*dt'[w(ic t t')-2Xm- 1 a s / 2 p^)f_ (55) 

The model's Eq.(53c) and the CQC Eq.(55) are identical if 

, 2Aa 3 / 2 2A S(x,i) 

W(x,t = S'(x,t) = - V y , 56a 

mo (pmorjlo 1 ' 2 W 

4(At)(^/o) = p(/3m c 2 ) 2 , (56b) 

where = £c 2 has ben used, and S'x,t) is the operator whose eigenvalues are s'(x, f) given in Eq.(54). 
Thus, this model produces CSL dynamics. 

It remains to check whether the numerical values in Eqs.(52), (56) are reasonable. Consider two possible temper- 
atures for the bath, the cosmic radiation temperature /3 _1 « 2 • 10~ 4 eV. and the Planck temperature /3 _1 s» 10 28 
eV. For normal mass densities, p w lgm/cc w 5 • 10 33 eV/cc, one gets respectively (3Gppa 2 w 10 -6 and 10 -38 , so the 
inequality of Eq.(52) is satisfied for both. 

For these two temperatures, (56b) gives p ~ 10 -112 and 10~ 49 respectively. The second number is too small for 
the conceptual picture we have outlined since, after 10 49 r « 10 5 sec, the spin in every cell is likely to have interacted 
and frozen, so the process would cease. However, one may change the conceptual picture. It is possible for a spin 
in a cell to interact repeatedly, without changing the dynamical equation because the bath states preserve the past 
orthogonality of the realizable states. Although they are represented by bath states, they are labeled by the past 
values of W(x,t). Each of these orthogonal realizable states at time t evolves independently of the other states over 
the next At, via the CSL evolution. 

Lastly, look at Eq.(56a). Note that W ~ S/N is intensive, proportional to the mean value of the spin in a cell, so this 
result is independent of the size of Ax, At as it should be. For the numerical coefficient in (56a), it is more informative 
to discuss S'(x, t)'s relation to S(x,t) rather than W(x,t)'s relation to either, because S'(x, t)'s mean value is simply 
p (see (53c)). One obtains from Eq.(56a), for a cosmic bath and for a Planck bath respectively, S' « (S'/A r )10 39 eV/cc 
and S' ~ (SyiV)10 71 eV/cc. Both factors are large compared to normal mass density ps 10 34 eV/cc, so that the mean 
spin value S/N for such a density is respectively 10 -5 and 10 -37 . 
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